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Abstract Transcranial Electrical Stimulation (TCES) and 
Deep Brain Stimulation (DBS) are two different applica- 
tions of electrical current to the brain used in different ar- 
eas of medicine. Both have a similar frequency dependence 
of their efficiency, with the most pronounced effects around 
lOOHz. We apply superthreshold electrical stimulation, spe- 
cifically depolarizing DC current, interrupted at different fre- 
quencies, to a simple model of a population of cortical neu- 
rons which uses phenomenological descriptions of neurons 
by Izhikevich and synaptic connections on a similar level 
of sophistication. With this model, we are able to reproduce 
the optimal desynchronization around lOOHz, as well as to 
predict the full frequency dependence of the efficiency of 
desynchronization, and thereby to give a possible explana- 
tion for the action mechanism of TCES. 

Keywords transcranial electrical stimulation, deep brain 
stimulation, integrate-and-fire model, desynchronization, 
high-frequency stimulation 



1 Introduction 

Deep Brain Stimulation (DBS) has attained much attention 
during the past fifteen years as a modem treatment of mis- 
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cellaneous neural and movement disorders, especially in the 
treatment of the symptoms associated with Parkinson's dis- 
ease (PD) 1 8 1, namely tremor, rigidity and bradykinesia, but 
also for epilepsy |39|, dystonia Illlll42ll23l . and essential tre- 
mor |7|. DBS also has promising effects in the treatment of 
obsessive-compulsive disorder (OCD) |31| , Tourette's syn- 
drom 017II4OI and depression ll34l . For the purpose of the 
treatment, electrodes are chronically implanted in a specific 
area of the brain. In the majority of cases of movement disor- 
ders, the targeted region is the subthalamic nucleus (STN), 
as the stalling of dopaminergic (inhibitory) neurons in the 
substantia nigra leads to pathological synchronized oscil- 
lations in the STN, which are coiTelated with the clinical 
symptoms of PD 1 13 1. Also other regions of the brain have 
been targeted for DBS, the globus pallidus internus (GPi) for 
the treatment of PD |32,6| as well as treatment of dystonia 
II42II23I . was attempted in the CAl region of the hippocam- 
pus 022II38I . and in the ventral intermediate nucleus of the 
thalamus (Vim) ||6|. Frequencies between 100 and 200 Hz 
(clinically often 130 Hz are applied) have been proven to 
give the best results in alleviating the symptoms, whereas 
low frequencies of roughly lOHz can even worsen the symp- 
toms. Understanding the working mechanisms of DBS from 
computational models is not an easy exercise: RecuiTent neu- 
ral networks with inhibitory neurons can exhibit lich behav- 
ior including synchronization 028114411411 : various mecha- 
nisms of desynchronization by stimulation lfT6l or nonlin- 
ear feedback |36|, to mention a few, have been proposed. 
Although there has been a lot of research, the action mecha- 
nism of DBS still remains elusive. 

At least in the last years, however, less effort has been 
made to analyze the action mechanisms of transcranial elec- 
trical stimulation, which is used for different purposes and 
is labelled with different names. Probably the most estab- 
lished one is Transcranial Electrical Stimulation (TCES), a 
technique which has been used to reduce drug requirements 
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for anaesthesia in surgical operations (TT]. Although it has 
been used in over 10000 operations at least up to 1998 [i27i . 
its working mechanism still remains to be explained. 

The TCES technique was developed by Limoge in the 
1970s f26l and the technical protocol as well as electrode 
placement is established under the terminus Limoge current 
lIZTl . Although the way of the Limoge current through the 
head is still not precisely known, and it is even not known if 
the current may, at least in parts, act by influencing periph- 
eral nerves outside the cranium [43], the part of the brain 
where transcranial electrical stimulation has its greatest ef- 
fect is presumably the neocortex, because current density 
decreases with the distance to the electrode. Since the neo- 
cortex is the part of the brain closest to the stimulation for 
all standard electrode positions, it is quite likely the part car- 
rying the greatest fraction of the current. Transcranial direct 
current stimulation in humans was also shown to enhance 
excitability f?",?! which on the one hand indicates plastic- 
ity effects and on the other hand shows that TCES can have 
a pronounced effect on neocortical regions and not only on 
specific subcortical structures as targeted by DBS. 

To account for an explanation of the action mechanism 
of TCES, we have built up a model of a population of neu- 
rons from the mammalian cortex, and applied external stim- 
ulation with different frequencies. This model is presented 
in section |2] In section [3] we give a short introduction to 
generalized phase definitions and introduce order parame- 
ters which can be used to quantify the phase synchroniza- 
tion of the neurons. Results are presented in section 2] and 
discussed in section |5] In the final section |6] we give a com- 
prehension and present possible extensions of the model. 



2 The model 

As we want to demonstrate the basic effects in a model that 
is as simple as possible, we do not attempt to account for 
biological parameter variability, as we expect that the basic 
mechanism does not rely on detailed properties of a specific 
cortical region. 



2.1 Neuron model 

Our model cortex consists of Nexc — 1024 excitatory and 
Ninh — 256 inhibitory neurons. This reflects the fact that 
the ratio of excitatory to inhibitory neurons (in the mam- 
malian cortex) is 4 : 1 191 1201 . Individual neurons were mod- 
elled according to the Izhikevich model of spiking neurons 
IfTSll . which contains two variables v and u, representing the 
membrane potential, and a recovery variable, e.g. a slow K+ 
current, respectively, v and u are governed by the differen- 



tial equations 

V = QMv'^ + 5v + 140 
ii — a{bv — u), 



(1) 
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where Text and /sj,„ represent the external and the synaptic 
input current, respectively. If u > 30, a reset is initiated, 



c, 



(3) 



Depending on the four parameters a, b, c, and d, the Izhike- 
vich neuron can map a rich variety of neuronal spiking pat- 
terns. We use a = 0.02, b = 0.2, c = -65, d = 8 for the 
excitatory subpopulation and a = 0.1, b = 0.2, c — —65, 
d = 2 for the inhibitory subpopulation. These parameter val- 
ues correspond to 'regular spiking' (RS) and the 'fast spik- 
ing' (FS) pattern, which are exhibited by most of the excita- 
tory and of the inhibitory neurons in the cortex, respectively 

Eina. 



2.2 Modeling synapses 

Each neuron makes a synaptic connection to any other neu- 
ron with probability p = 200/(A*'e2:c + Ninh), so each neu- 
ron has 200 synaptic connections on average. In reality, there 
are thousands of synapses per neuron [20] , whereof here 
we model a small subnetwork of a local cortical assembly 
and therefore cannot explicitely consider long-range con- 
nections. Consequently, we also neglect structured connec- 
tivity, apart from taking into account a sparse random con- 
nectivity. The synaptic input current to a neuron i is given 
by the sum over the postsynaptic potentials of all neurons j 
presynaptic to i 
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where Xj and Xi are the types (i.e. excitatory or inhibitory) 
of the pre- and the postsynaptic neurons, respectively, and 
can take on the values E and /. Here (ij,i) is the adja- 
cency matrix describing whether synaptic connections be- 
tween j and i exist. The sign variables s^ formally account 
for the excitatory or inhibitory nature of the presynaptic neu- 
ron, with se — +1 and sj = —1. Then, the gx ,xi are 
the synaptic strengths between two types of neurons, with 
gE,E ^ 0.6, QEj = 0.1, gi^E = 0.2 and gij ^ 0.05. We 
have adapted the ratios of these values according to the ref- 
erence model [101 . For the axonal delay we use 5 — 0.25ms; 
Tj; is a time constant depending on the presynaptic neuron, 
with 0.2ms for excitatory and 0.4ms for inhibitory neurons, 
respectively, and tj denotes the time since the last spike of 
the presynaptic neuron j (in ms). 

To consider the dynamics of a local neural assembly un- 
der stimulation, it is of low relevance from where incom- 
ing projections originate. To ensure non-trivial activity in 



Desynchronizing effect of high-frequency stimulation in a generic cortical network model 



the network, external input to the neural network is assumed 
from Nfixt = 128 external excitatory neurons. These neu- 
rons are not explicitly modelled, but they make synaptic 
connections into the network with the same properties as 
the network synapses, and fire action potentials at random 
with probability 0.01 during each time step. The model was 
integrated using a fourth-order explicit Runge-Kutta method 
with a time step of 0.05ms. 



3 Phase synchronization 

To quantify synchronization of an ensemble of oscillators, 
more convenient observables than correlation functions can 
be used. Here we follow an established approach 1171 1241 to 
define a complex-valued order parameter based on a gener- 
alized phase in a phase space appropriately chosen for the 
model at hand. 



3.2 Order parameter 

We define the complex- valued order parameter f of the phases 

as 
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This definition is according to the definition of the order pa- 
rameter in the Kuramoto model of phase synchronization |[T] 
l24il . We analyze the synchronization in the whole population 
(fw), the excitatory subpopulation (f^) and the inhibitory 
subpopulation (f/) as well, so N can be N^xc + Ninh, N^xc 
and Ninh, respectively. The absolute value r of f lies in the 
interval [0, 1], where corresponds to a complete unsyn- 
chronized state and 1 to complete synchronization. 

As these order parameters oscillate already for the un- 
stimulated system, we are mainly interested in the average 
of the order parameters, which we denote by a bar, e.g. rw- 



3.1 Generalized phase 

There are at least three concepts for a generalized phase 
definition of chaotic oscillations fFS\. For the sake of nu- 
merical simplicity, we use a geometrical phase angle defi- 
nition. In absence of external stimulation, each neuron has 
an attracting limit cycle in the {v, u, v, u}-space; to define 
a phase, this trajectory has to be suitably projected into a 
two-dimensional plane, which indeed is possible here. We 
define the phase of the oscillator i (i.e., the neuron i) as the 
angle between a given (fixed) direction and the position of 
the neuron's state in a chosen plane of the phase space. 



i{t) — arctan 
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where the point (wi,c, —Vi,c) is within the rotation centre of 
the attractor of the neuron i in the plane of the phase space 
spanned by Vi and ~Vi, with Vi being the time derivative of 
the membrane potential Vi of the neuron i. 

The first choice for such a plane in the phase space to 
be used for phase definition would be the {v, u}-space, but 
unfortunately the attractor in this plane has multiple rotation 
centres and changes its position for different stimulation in- 
tensities, so the {v, — v}-space is a better choice. To ensure 
a monotoneous increase of 0, one has to choose a coordinate 
centre within the trajectory loop; here we select the point 



{Vi. 



~) = (ci, 0) as rotation centre. We further choose 



{v, — v} (rather than {v, v} which results in </) < 0) to en- 
sure > 0, for convenience. Other linear combinations of 
V, u and their first time derivatives could be used as well as 
an embedding of the dynamics. 



4 Results 

4. 1 Behaviour of the unstimulated system 

In the unstimulated system, we observe tightly synchronized 
spike volleys that occur with ^ lOHz. This is in the alpha- 
range which is associated with an alert but relaxed state. 
Iil4i . This synchronization phenomenon is widely known for 
random coupled neural networks [2, 18. 21, , 91. 

As one can already suspect by looking at the distribution 
of the spikes over time, all three order parameters have big 
values close to 1, at least between spike volleys. When such 
a population spike occurs, the order parameters break down 
and reach values that can be as small as 0.2 or even smaller 
for a very short time. 

The dynamics of the unstimulated system is shown in 
Fig.[T]in each of the subfigures (a-h) before the stimulus on- 
set. At 3000ms, stimulation impulses were applied accord- 
ing to different protocols as described below. For each stim- 
ulation type. Fig. |2] shows the respective network activities 
(left panels), as well as the time-dependence of the order 
parameter (right panels). 



4.2 Uninterrupted DC stimulation 

For (uninterrupted) DC stimulation, the frequency of the 
spike volleys is increased, but the spike volleys of the excita- 
tory subpopulation still remain discriminable. In contrast to 
that, the inhibitory subpopulation now fires continuously, as 
there are only few and relatively weak inhibitory-inhibitory 
synapses which could affect this subpopulation to cease fire. 
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Fig. 1 Spike histograms of 3000ms of four simulations, for the excitatory (left-hand side) and the inhibitory (right-hand side) subpopulation. The stimulation 
protocol is as follows: Stimulation lext was switched on 3000ms after initialization for 1000ms. Each stimulation is monopolar (excitatory DC stimulation) 
and either uninterapted DC (a,b), or modulated by a rectangular envelope of lOHz (c,d), lOOHz (e,f), or lOOOHz (g,h), respectively. Before and after the 
stimulation, the behavior of the unstimulated system is visible. 
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Fig. 2 Left-hand side: Spike times of the neurons in the network. Each point denotes a spike. The ordinate denotes the number of the neuron if numbered 
serially. Above the black line are the inhibitory neurons, below it the excitatory neurons. Right-hand side: The order parameter tiy over the time (ordinata: 
relative to stimulation onset t = 0). Amplitude of stimulation is 10 in both cases. Note the different scaling of the abscissae on the left- and the right- 
hand side. (a),(b) before stimulation. (c),(d) uninterrupted stimulation. (e),(f) lOHz rectangular pulses. (g),(h) lOOHz rectangular pulses. (i),(j) lOOOHz 
rectangular pulses. 
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For this type of stimulation, Tw is decreased. In fact, it 
is the strongest reduction ofrw that we observe for all stim- 
ulation frequencies. But if one takes a view on the order pa- 
rameters of the two subpopulations, one gets a different find- 
ing. Te is even increased. In contrast, rj fluctuates strongly 
and fast almost between 1 and 0. Therefore, fj is decreased 
all in all, but not as strong as rw- So, the reduction ofrw 
stems in parts from the desynchronization of the inhibitory 
subpopulation, but also from the desynchronization of both 
subpopulations from each other. 



4.3 Low frequency stimulation (lOHz) 

For low frequency stimulation, the spiking pattern of the 
excitatory subpopulation remains almost unchanged, but as 
such frequencies are near the systems eigenfrequency, we 
have a resonance phenomenon. Therefore, the population 
spikes are now time-locked to the stimuli, and occur with 
exactly the freqency of the stimulation. Hence the activity 
becomes more synchronized, leading to the observable ef- 
fect that all three order parameters are increased (Fig. |33. 



4.4 High frequency stimulation (lOOHz) 

If we increase the frequency of the stimulation to high fre- 
quencies of ^ lOOHz, the spike volleys still occur during 
stimulation and synchronize with the stimulation frequency. 
If the amplitude of the stimulation is not high enough, some 
stimuli may be missed, so that, for example, three spikes oc- 
cur locked to the stimuli, and then for another two stimuli, 
the neural network is silent. 

This finding is supported by the order parameter rw, 
which decreases around lOOHz. It does not decrease as much 
as for the uninterrupted stimulation, but in contrast to that 
kind of stimulation, the order parameter te of the excita- 
tory subpopulation decreases as well, and below the baseline 
value of the unstimulated case. Similarly, r/ is decreased, 
but the minimum is not as sharp as for Te- 



4.5 Very high frequency stimulation (lOOOHz) 

If we increase the stimulation frequency further, the indi- 
vidual neurons cannot follow the stimulation frequency, and 
we observe spike time characteristics that are very similar to 
the unstimulated case. All three order parameters again in- 
crease to higher values. One could think that it is possible to 
reproduce the desynchronizing effect of the moderate high 
frequencies by increasing the stimulus amplitude, but unfor- 
tunately we could not generate such effect (not shown). 



4.6 Post-stimulation characteristics 

As we did not model aspects like synaptic plasticity, there 
are no long-term effects of the stimulation. All effects cancel 
out with cessation of stimulation, except a small reduction in 
spike rate, which soon returns to baseline values. 



5 Discussion 

We applied depolarizing electrical stimulation to a simple 
cortex model. The stimulating current is interrupted with 
different frequencies. 

If not stimulated, the model exhibits synchronized spike 
volleys or population spikes at frequencies that lie in the 
alpha range. This synchronization phenomenon is widely 
known for randomly coupled neural networks (|2]|9][T8]|2T] . 
The three order parameters we have defined are close to the 
maximum value 1 when there are no spike volleys, but break 
down when population spikes occur This is not surprising as 
the movement on the trajectories is very fast during spikes, 
compared to the time between spikes. Thus, if there is only 
a little difference between the spike timings of two different 
neurons, they may be far apart from each other, whereas for 
the time between population spikes, all neurons have values 
of V and —v located in a small volume of the phase space. 
Thus the phases are very similar between population spikes 
and more different during the population spikes, and the or- 
der parameters break down in the latter case. 

For uninterrupted DC stimulation the frequency of the 
spike volleys increases, and the inhibitory subpopulation be- 
gins to fire continously. Although the average order parame- 
ter rw of the whole population decreases strongly, the order 
parameter Te of the excitatory subpopulation is increased 
by uninterrupted stimulation. Enhanced synchronization in 
large groups of neurons is also reason for epileptic spasms. 
As motor neurons (which are - excitatory - pyramidal neu- 
rons) can be effected by TCES as well as any other popu- 
lation in the cortex, we think this effect could be a possible 
explanation of the clonic spasms which are produced by un- 
interrupted high intensity extracranial DC stimulation [,26l . 

If the stimulation is interrupted at low frequencies, pop- 
ulation spikes become time-locked to the stimuli, and the 
order parameters Tw and Te even increase. This is not un- 
expected as stimulation near the eigenfrequency of a sys- 
tem leads to resonance. If we think about DBS, for which a 
similar frequency dependence of effectiveness as for TCES 
is observed (e.g., compare |I26]|27]|33 for TCES with ^ 
[30l for DBS), it is known that stimulation with low frequen- 
cies can even worsen the symptoms which are associated 
with too strong pathological synchronization lfT2l . In terms 
of TCES, such a stimulation would possibly counteract any 
anaesthetic action, but this cannot be said definitely as there 
seem to be no experiments on that yet. 
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Fig. 3 Dependence of the order parameter on the frequency of stimu- 
lation. Intensity of stimulation is 10. (a) ryy of the whole population, 
(b) r_B of the excitatory subpopulation. (c) r/ of the inhibitory subpop- 
ulation. In all three plots, averages are shown over 15 simulations, and 
the corresponding standard deviations are indicated by error bars. 



For stimulation with moderately high frequencies, we 
found the introduction of a new activity pattern. As the neu- 
rons try to follow the stimulation frequency, smaller spike 
volleys occur which are locked to the stimuh. But due to the 
short duration of the stimuli, some of them are missed by 
the population spikes, and only a smaller number of neurons 
participates in each such spike volley, compared to the in- 
trinsic volleys. Effectively, this is a desynchronizaton of the 
neural activity. It is well known that Deep Brain Stimula- 
tion of the subthalamic nucleus with such moderately high 
frequencies cancels out the pathological synchronization in- 
duced by Parkinson's disease. We propose that in the cortex 
this desynchronization phenomenon intercepts neural infor- 
mation processing, as it is known from experiments on elec- 
troanaesthesia (EA) Il26l . which eventually lead to the in- 
vention of TCES. 

For very high stimulation frequencies, the desynchroniz- 
ing effect again vanishes. This is in compliance with investi- 
gations in DBS |22J, reporting stimulation in the kHz range 
to have no effect, and with old reports on EA showing the 
current intensity necessary for the loss of consciousness to 
greatly increase for such high frequencies [35 1. In fact, in- 
creasing the stimulation amplitude could not fully compen- 
sate the effect of increasing the frequency in our model. In- 
stead, by increasing the amplitude of stimulation, we rather 
found effects similar to those of uninterrupted current. Per- 
haps this explains why currents with a frequency of approxi- 
mately lOOHz, modulated with ultra high frequency currents 
{r^ lOOkHz), are comparatively effective than the unmodu- 
lated current 026II27I . If the ultra high frequency current acts 
similar to an uninterrupted current, the effect is expected 
to be similar to the unmodulated lOOHz current, but it has 
lesser secondary effects, as less current is delivered to the 
brain. 

6 Summary and Outlook 

Deep Brain Stimulation and Transcranial Electrical Stimu- 
lation are two promising applications of electrical currents 
in medical practice. Whereas DBS has largely replaced pal- 
lidotomy in the treatment of Parkinson's disease [i30l . and 
is prosperous in the treatment of some other neural diseases 
ir7l fm[3n[39ll42l . TCES yet does not have reached a com- 
parative level of potential applications that it could have in 
surgical practice. 

In summary, we have presented a generic cortex model 
which reproduces the frequency dependence of the activity 
of TCES: With uninterrupted stimulation, the overall syn- 
chronization is decreased but the synchronization of the ex- 
citatory subpopulation is increased. For moderately high fre- 
quencies synchronization of the whole population and the 
excitatory subpopulation is decreased as well. We assume 
the same mechanism to account for the similar frequency 
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dependence of the efficacy of DBS, resulting in a drastic 
desynchronization of the spiking activity. Stimulation of our 
cortex models with unipolar lOOHz rectangular waveforms 
leads to distinct changes of the spiking properties of the neu- 
ral network, and it markedly reduces the synchronization 
of the network, which is quantified by the order parame- 
ters rw, te and r/. Whereas slower or higher stimulation 
frequencies do not lead to such distinct changes, and in- 
stead can even tend to enforce the natural behaviour of the 
model, moderately high frequencies (around 100 Hz) lead to 
a desynchronization of activity. We therefore suppose that 
desynchronization of cortical activity and the introduction 
of cortical noise is at least in parts accountable for the ef- 
fects of TCES, as it disrupts ongoing signal processing in 
the cortex. 

This model obviously just gives hints to the the question 
how TCES actually works. To consolidate our theory that 
the desynchronization of cortical activity leads to the ben- 
eficial effects of TCES we would have to build up a more 
complex and of course larger model containing more neu- 
rons, which should take into account the spatial structure 
of the brain. If a spatial distribution of the neurons is mod- 
eled, one should desist from assuming homogenous electri- 
cal fields and model explicitly the electric field distribution. 
Also, neurons could be modeled in different degree of detail, 
e.g., according to Izhikevich ifTSl taking into account more 
firing patterns, or as advanced versions of Hodgkin-Huxley 
neurons with realistic synaptic dynamics, including synap- 
tic plasticity and axonal delays that depend on the spatial 
distances of the neurons. Similarily, a multi-compartment 
model - containing at least two or three compartments for 
axon, soma, and eventually dendrite - could be used. While 
the quantitative dynamics of such more detailed models might 
come closer to reality, we expect those models to exhibit 
a similar desynchronization effect in the lOOHz regime, as 
demonstrated in our model. 
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